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ABSTRACT 


Some  exact  solutions  for  the  time-dependent  incompressible 
irrotational  flow  of  a fluid  contained  between  two  free  surfaces 
are  presented.  The  free  surfaces  are  either  concentric  circular 
cylinders  or  concentric  spheres.  The  cases  where  constant  pres- 
sures are  applied  at  one  or  the  other  of  the  free  surfaces  are  also 
described,  and  the  corresponding  flows  are  analyzed  for  Taylor 
instability  of  the  free  boundaries.  When  the  free  boundaries  are 
concentric  circular  cylinders,  flows  with  constant  non-zero  circu- 
lation for  circuits  containing  the  cylinder  axis  are  also  described. 
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1.  INTRODUCTION 


The  exact  hydrodynamic  flows  described  in  this  report  were  found 
incident  to  carrying  out  the  analytical  phase  of  a project  entitled 
"Ship-Wave  Interactions"  funded  by  the  Office  of  Naval  Research 
under  Task  No.  NR  334-003.  So  far  as  we  know,  these  exact  solutions 
have  not  been  reported  elsewhere.  We  will  say  a little  about  our  moti- 
vation for  finding  these  solutions  in  the  next  paragraph,  but  our  reasons 
for  presenting  them  here  in  a separate  report  can  be  stated  very  simply: 
To  put  before  the  numerical  hydrodynamics  community  some  simple,  but  non- 
trivial, exact  solutions  of  the  non-linear  equations  for  the  motion  of 
a liquid  with  a free  surface,  which  might  be  used  to  check  the  accuracy 
of  various  numerical  techniques  for  solving  time-dependent  hydrodynamic 
free  boundary  problems. 

Our  reason  for  introducing  the  flows  described  here  is  an  outgrowth 
of  a relation  that  has  been  found  between  turbulence  and  energy  non- 
conservation for  inviscid  flows.  Details  of  that  connection  will  be 
reported  elsewhere  (Ref.  1).  Briefly,  for  a hydrodynamic  flow  it  turns 
out  that  a convenient  measure  of  the  local  turbulence  of  the  flow  is 
given  by  the  quantity 


Pou3)  + 7 • pou^p)]  , 


(1.1) 


where  u(x,t)  is  the  velocity  field,  p0  is  the  constant  liquid  density,  and 
P is  the  pressure.  The  presence  or  absence  of  turbulence  is  characterized 
by  whether  or  not  Q < 0 (more  precisely,  in  a stochastic  framework,  by 
whether  or  not  Q < 0 "almost  surely"). 
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In  the  next  paragraph,  we  will  give  an  example  of  a flow  with  u • i ® 
u • j = 0 and  u * k a. function  of  z and  t only,  for  which 

Q oc  - 6(z-z<,)  6(t-to)  . (1.2a) 

A natural  question  which  arises  in  this:  Is  there  a hydrodynamic  flow  with 

— ♦ — ♦ -4 

u * k = 0 and  u * i,  u • j independent  of  z,  such  that 

Q « - 5 (x-Xq)  6(y-y0)  6(t-to)  , (1.2b) 

and  is  there  a flow  with 

Q “ “ 6(x-Xo)  6(y->'c)  6(z-Zo)  f (t-to)  ? (1.2c) 

Further,  we  may  ask  if  somehow  one  car.  exore'a  gerarcl  turbulent  flows  for 
different  numbers  of  independent  spec ? a:  «!>'*«  in  terms  of  superpositions 
of  non-turbulent  flows  and  eleme*  *•* .:  - * .ows  with  turbulence  characterized  by 
(1.2a),  (1.2b),  or  (1.2c).  At  this  -vicing,  we  have  not  answered  the  latter 
inquiry,  and  indeed  there  will  be  no  need  to  answer  it  until  our  analysis 
has  proceeded  to  a more  advanced  stage.  As  regards  the  other  question,  in 
this  report  we  construct  flows  having  properties  (1.2b)  and  (1.2c)  as 
natural  generalizations  of  the  flow  with  Q satisfying  (1.2a). 

For  a flow  satisfying  (1.2a),  consider  the  hydrodynamic  flow  with  ini- 
tial conditions 

/ 0 0 £ | z|  < Zq 

p(x,0)  = / p0  Zq  £.  |z|  £ zx  (1.3a) 

( 0 zx  < | z| 
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u (x,0)  - - U k agn(z)  , z0  S |z|  £ zx  . 


(1.3b) 


Here  0>  0.  For  0 < t < we  have 


P(x,t) 


0 0 £ |z|<  (zq  - Ut) 

Po  (zo  - Ut)£  |z|  £(z,,  - Ut) 

0 (zx  - Ut)<  |z| 


(1.4a) 


u (x,t)  = - U k agn(z),  (Zq  - ut)  S | z|  s (^  - Ut) 


(1.4b) 


For  later  times,  we  choose  the  flow  that  satisfies 


P(x,t) 


Po  0 £ 1 z|  £ (zx  - Zq) 

0 (Zi  - Zo><  |z| 


(1.5a) 


I 

l 

I 

I 

i 

i 

i 

II 


u (x,t)  - 0,  0 £ |z|  s(z1  - Zq).  t > 


d.5b) 


For  the  flow  given  by  (1.3) -(1.5),  we  readily  find 


Q - - Po  ^(zi-zo)  5(z)  6(t-Tf) 


(1.6) 


This  is  a flow  of  type  (1.2a);  others  may  be  found  by  shifting  the  origins 
of  space  and  time. 

Let  us  make  the  following  observation:  Instead  of  (1.5),  one  might 


choose  the  flow 


p(x»t) 


D 0 £ | z | < U(t-  -j^) 

Po  U(t-  -22.)  £ |z|  s;  z^  Zq  + U(t-  ■^•),  t * ^ , (1.7a) 

0 zx  - zq  + U(t-  ^)  < | z| 
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and 

u (x,t)  - U k sgn(z),  U(t-  “)  * U1  ^ z1-  Zq  + U(t-  -y)  , t > . (1.7b) 

For  the  flow  given  by  (1.3),  (1.4),  and  (1.7),  Q = 0.  We  may  think  of 
the  flow  (1.3)-(1.5)  as  "perfectly  inelastic",  the  flow  (1.3),  (1.4),  (1.7) 
as  "perfectly  elastic".  There  are,  of  course,  various  intermediate  possi- 
bilities. The  reason  we  have  focused  on  the  flow  (1.3)-(1.5)  is  that  this 
ia  the  flow  that  will  be  calculated  by  our  water-wave  algorithm  (Ref.  2). 

That  is  not  a mandatory  feature  of  our  calculation,  but  it  seems  like  a 
good  first  approximation,  and  it  is  with  that  feature  in  mind  that  we  have 
analyzed  stability  and  turbulence  for  the  flows  generated  by  the  algorithm. 
However,  other  numerical  schemes  may  calculate  instead  the  flow  (1.3),  (1.4), 
(1.7),  or  possibly  intermediate  cases.  The  exact  flows  we  present  in  this 
report  are  not  predicated  on  the  assumption  of  inelasticity  of  "collisions" 
and  will  be  of  use  whatever  numerical  scheme  one  may  use.  The  only  difference 
will  be  in  how  the  various  flows  presented  here  are  to  be  "patched  together" 
at  moments  of  "collision". 

Like  the  flow  (1.3)-(1.5),  the  flows  to  be  presented  in  the  next  sec- 
tion have  a free  boundary  consisting  of  two  components  unconnected  to 
each  other.  Furthermore,  they  are  "one-dimensional"  in  the  sense  that  the 
velocities  depend  on  only  one  spatial  coordinate,  if  one  uses  cylindrical 
or  spherical  coordinates.  In  most  of  the  cases,  there  is  only  one  non-zero 
velocity  component  in  an  appropriate  coordinate  system.  All  the  flows  are 
irrotatlonal  as  well  as  incompressible.  As  a consequence,  it  turns  out  that 
the  "cylindrical"  solutions  have  a greater  variety  than  the  "spherical"  ones 
because  of  the  possibility  of  having  non-zero  circulation  for  an  irrotatlonal 
flow  in  a multiply-connected,  but  not  simply-connected,  domain. 
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There  are  no  rigid  boundaries  for  the  flows  presented  here,  and  in 
addition  the  domains  occupied  by  the  moving  liquid  are  bounded.  It  has 
been  seen  (Ref.  3)  that  under  such  conditions  the  effect  of  gravity  is 
only  to  bring  about  an  acceleration  of  the  whole  flow  configuration  in 
the  direction  of  the  gravitational  attraction.  Thus,  there  is  no  loss  of 
generality  in  ignoring  the  effect  of  gravity  in  these  flows.  If  there  is 
gravity  present;  the  flows  should  be  considered  to  be  given  in  an  appro- 
priately accelerated  system  of  coordinates.  Of  course,  if  one  desires  a 
more  stringent  test  of  his  solution  algorithms,  not  only  may  he  compute 
the  flow  in  the  presence  of  gravity,  but  he  may  solve  the  equations  in 
a coordinate  system  for  which  the  special  solutions  given  here  do  not 
possess  such  a simple  symmetry.  In  addition  to  accelerated  systems, 
such  coordinates  may  be  obtained  by  displacing  the  origin  of  coordinates 
and  moving  the  system  in  uniform  rectilinear  motion.  Since,  however,  the 
solutions  given  here  are  all  relatively  smooth,  the  demands  made  on  a 
solution  algorithm  to  duplicate  them  will  be  correspondingly  limited, 
and  the  test  on  the  algorithm  afforded  by  comparison  of  the  exact  solu- 
tion with  the  approximate  one  will  be  one  of  necessity,  not  of  sufficiency. 
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2.  DESCRIPTION  OF  THE  FLOWS 

For  all  the  flows  given  here,  the  liquid  occupies  the  domain 

. ro  (t)  < r < n(t)  , (2.1) 

where  for  three-dimensional  flows  r = (x2 +y2+z 2)^  and  for  two-dimensional 
flows  r = (3^+y2)^.  The  sets  {x  | |x|  = r0(t)}  and  [x  | |x|  = rx(t)}  are 
called  the  inner  and  outer  free  surface,  respectively.  The  cases  we  con- 
sider are;  (i)  three-dimensional  flow,  zero  pressure  on  the  free  surfaces; 

(ii)  two-dimensional  flow,  zero  pressure  on  the  free  surfaces;  (iii)  three- 

dimensional  flow,  pressure  applied  at  one  of  the  free  surfaces;  (iv)  two- 

dimensional  flow,  pressure  applied  at  one  of  the  free  surfaces;  (v)  two- 

dimensional  flow  with  circulation,  zero  pressure  on  the  free  surfaces;  and 
(vi)  two-dimensional  flow  with  circulation,  pressure  applied  at  one  of 
the  free  surfaces. 

We  caution  the  reader  that,  for  economy  of  notation,  the  same  symbol 
may  have  different  meanings  in  the  descriptions  of  different  flows.  Of 
course,  if  an  equation  appearing  in  one  flow  description  is  referred  to 
in  another  flow  description,  the  symbols  that  appear  therein  will  be  con- 
sistent. 

THREE-DIMENSIONAL  FLOW,  ZERO  PRESSURE  ON  THtf  FREE  SURFACES 

Since  the  flow  is  incompressible  and  irrotational,  the  velocity  field 
is  derivable  from  a potential  cp  satisfying  A<P  * 0.  Bernoulli's  equation 
is 

<Pt+lW2+~  = 0 , (2.2) 


where  P is  the  pressure  and  po  is  the  constant  liquid  density. 


: ! i*h 


s . 


7 


hi 


.. 

t i 


0 


i n 
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We  write 


cp  = A(t) 


(do-T) 


+ B(t) 


(2.3) 


There  are  four  conditions  to  determine  A(t),  B(t),  r0(t),  and  rx (t) : The 
two  kinematical  conditions  that  points  on  the  free  surfaces  remain  there, 


• A , • A 

ro  = and  rx  = -jj  , 


(2.  A) 


and  the  two  dynamical  conditions  that  the  pressure  on  the  inner  and  outer 
free  surfaces  be  zero.  Using  these  pressure  conditions  we  get,  on  substi- 
tuting equation  (2.3)  into  (2.2), 


•A  1\  A • ■ , 1 A2 

A VU  *?'•  +B+i^-° 


(2.5) 


and 


A • • 1 A2 

- ~3  ri  B + 2 t = 0 

n ri 


(2.6) 


From  the  kinematical  conditions  (2.4),  we  get  conservation  of  volume: 


r®(t)  - r£(t)  = t?(0)  - rg(0)  = V . 


(2.7) 


By  combining  equations  (2.4)-(2.6)  we  also  get  conservation  of  energy: 


*a<t>  (mo  • mo)  ■ k‘m  Odor  • * k)  ■ E • (2-8) 


It  is  easier  to  work  with  the  equations  in  dimensionless  form.  In  terms 
of  E and  V above  we  can  find  a characteristic  time 


3 

I 


1 


t 
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* v^6 

c =?7r*- 


(2.9) 


Letting 


t n ro 

T = 7*  » - ^77 


*1  = -77T  » and 
IT* 


gi/s  yi /e 


(2.10) 


we  get 


R?  - R^ 


(2.11a) 


_1 1_ 

Ro  Rj. 


c*2  ’ 


(2.11b) 


where 


(2.11c) 


(2. lid) 


If  we  use  (2. lla,b,d)  to  express  a in  terms  of  Vx  , we  see  that  (2.11c) 
is  an  equation  for  t in  terms  of  V.  Let  us  consider  a special  solution  for 
which,  at  the  origin  of  time,  t = 0,  we  have  Vx (0)  = 1,  or  Rq (0)  = 0.  Then 
we  have 


T.i 

1 3 

J 1 


(S-i)v* 


(2.12) 


There  are  two  possible  signs  for  the  square  roots  in  (2.12).  We  shall 
specialize  to  the  case  where  the  positive  signs  are  taken. 
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We  note  the  following  limiting  values  of  the  integral  in  (2.12). 
When  Vx  1, 

T - | -1)B/6  . 


(2.13) 


Equivalently, 


Vl 


„ .6/6 

(*T) 


+ 1 , 


(2.14a) 


Ro  « (Vx-l)l/3  - ( | r'fB 


(2.14b) 


On  the  other  hand,  when  Vx  “*  00 , 


1 l/a 

T - — Vx  , 

n/3 


(2.15a) 


Rx  -•  \/3  T 


(2.15b) 


In  dimensional  variables, 


ro 


*/B(i 


. s/e 


0 


as  t I 0 


and 


r?  E\l/s 

^ Iv!  *=  as 


t t 


(2.16a) 


(2.16b) 


Apart  from  limiting  values,  we  see  that  (2.12)  gives  R:  as  a monotonically 
increasing  function  of  T,  which  we  denote  by 
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Now,  suppose  we  have  a general  Initial-value  problem,  with  ro(0), 
ri (0) , and  ri<0)  given.  From  (2.4),  (2.7),  and  (2.8)  we  compute  V and 
E.  With  Rx  (0)  obtained  from  rx  (0)  and  (2.10),  we  use  (2.17)  to  compute 
a ''time"  T0  > 0 such  that 


MO)  - F(t0)  . 
If  rx (0)  > 0,  the  flow  for  t > 0 is  given  by 


(2.18) 


rl  (t)  - V 1/3 


(2.19) 


if 

where  t is  given  in  (2.9).  r0(t)  is  found  from  (2.7).  If  rx (0)  < 0,  the 

flow  for  0 < t < T0t*  is  given  by 


ri  (t)  =W3  F (r0  - fff)  , 


(2.20) 


with  r0(t)  given  again  by  (2.7).  For  t > T0t*,  the  case  of  "inelastic 
collision"  gives 


ri  (t)  = V173 


(2.21a) 


and  in  the  case  of  "elastic  collision" 

rl(t)  -V1*  F - To)  . 


(2.21b) 


In  the  former  case,  with  Q given  by  (1.1),  we  have 


Q - - 2tt  oo  E 6(t-T0t* ) 6(x)  . 


(2.22) 


In  the  latter  case,  Q = 0. 


- — W 


MNllMUa 
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TWO-DIMENSIONAL  FLOW,  ZERO  PRESSURE  ON  THE  FREE  SURFACES 

The  similarities  with  the  three-dimensional  case  are  so  extensive  that 
we  only  indicate  what  changes  need  to  be  made  in  the  equations  above  as  we 
go  along. 

In  place  of  (2.3)  we  have 

<p  = A(t)  in  + B(t)  . (2.23) 

The  kinematical  boundary  conditions  are 

r0  = ~ and  rx  « ■£-  , (2.24) 

rO  rl 

from  which  we  get  conservation  of  volume, 

r?  (t)  - r§(t)  = t?(0)  - r|(0)  - V . (2.25) 

The  conservation  of  energy  is  expressed  by 

2”§7iy-AS<°>  • <2-26> 


The  characteristic  time  is 
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Then  the  equations  assume  the  dimensionless  form 


where 


R?  - Ro  - l , 

(2.29a) 

A R1  _ 1 

11  Ro  " 02  * 

(2.29b) 

dy., 

= 2“  ’ 

(2.29c) 

Vi=  R?  . 

(2.29d) 

Consider  the  flow  for  which  V* (0)  = 1,  Rq(0)  - 0.  In  place  of  (2.12) 
we  get 


T 


(2.30) 


where,  as  above,  we  take  the  positive  sign  for  the  square  root.  This  gives 
Rx  as  a monotonically  increasing  function  of  t,  denoted  by 


Rx  * F(t)  . 

The  following  limits  are  of  interest.  When  Vx  -*  1, 


(2.31) 


(2.32) 


(2.33a) 


I 

I 

1 

I 

1 
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When  Vj  "*  ® , 


1 l /a 
T-~V1 


Ri  -*  T . 


In  dimensional  variables, 


23/a  gi  A ti/a 


In 


as  t i 0 


and 


p < 

"*  **  ( v ) c as  t t “ 


15. 


(2.34a) 


(2.34b) 


(2.35a) 


(2.35b) 


Suppose,  in  general,  we  are  given  ro(0),  rx(0),  and  ^(0).  From  (2.24)- 
(2.27)  we  compute  V,  E,  and  t*.  Find  Rx  (0)  by  (2.28),  and  then  use  (2.31) 
to  find  t0  > 0 such  that 


Ri(0)  * F(t0)  • 


(2.36) 


For  rx  (0)  > 0,  the  flow  for  t > 0 is 


MO  - V1/3  f(^0+  pj)  . 


(2.37) 


For  rj (0)  < 0,  the  flow  for  0 < t < t*  t0  is 


• • 


MO  ’ F (o-pf)  ‘ 


(2.38) 
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For  t > t*  t0  and  r.  (0)  < 0,  we  get 


ra(t)  * W2 

In  the  inelastic  case,  and 

ri  (t)  - W2F^Lf  - t0^) 

in  the  elastic  case.  In  the  former  case. 


(2.39a) 


(2.39b) 


Q = - tt  p0  E 6(t-t*T0)  6 (x)  . (2.40) 

THREE-DIMENSIONAL  FLOW,  PRESSURE  APPLIED  AT  ONE  OF  THE  FREE  SURFACES 

The  velocity  potential  is  given  by  (2.3),  and  the  kinematical  conditions 
at  the  free  surface  are  still  given  by  (2.4).  Since  the  only  quantity  of 
interest  is  the  difference  between  the  pressures  applied  at  the  inner  and 
outer  surfaces,  we  have  no  loss  of  generality  in  setting  the  pressure  at 
the  inner  surface  equal  to  zero.  Let  the  pressure  at  the  outer  surface  have 
the  value  Px  constant  on  the  surface.  Then  equation  (2.5)  scill  holds,  but 
(2.6)  should  be  replaced  by 


A.-1A2*!  _ 

W'1  +B+2^  + ^’° 


(2.41) 


The  conservation  of  volume  is  expressed  by  (2.7).  The  conservation  of  energy 
takes  the  following  form  when  Px  is  Independent  of  the  time: 


A2<t)(r°(t>  " ri  (t))  + 3 p0  r°(t)  “ A2(0)(ro(0)  " r1(0))+  3 p0  r°(0)  = 6 


(2.42) 
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Consider  the  case  when  e,  as  defined  by  (2.42),  is  > 0.  This  is 
the  case  whenever  Pj^  > 0 and  sometimes  when  Px  < 0.  Then  let  us  write 
e = E.  Define  t*,  t,  Rq,  Rx  , or,  and  Vi  by  (2.9),  (2.10),  and  (2. lid). 
Define  TTX  by 


V h 
E p0 


(2.43) 


Then  the  governing  equations  are  (2.11a),  (2.11c),  and 


(2.44) 


(2.45) 


Ue  take  the  positive  signs  for  the  square  roots.  Hie  behavior  for  small 
time  (t  1 0)  is  given  by  (2.14a,b). 

Two  cases  should  be  distinguished.  In  the  first,  Px  > 0.  In  that 
case  the  expression  (2.45)  for  t fails  for 


Vi  > §o 


(2.46) 


- A-J.n. 
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Define  t*  by 


i y ('-mtt 

1 (£-l)1/e  (l-  f nx' (S-l)J 


_*  l 


Then  as  (t*-t)  l 0,  we  get 


v>  - 1 + 2^r  - ! "> 


3 Kr  &r 

- c»*r 

\0a(^r<0i 


(T*-T)2 


(T*  - T)2 


(2.47) 


(2.48a) 


(2.48b) 


For  0 < t < t*»  (2.45)  gives  Rx  as  a monotonically  increasing  function  of  t, 


which  we  denote  by 


Ri  - F(t) 


(2.49) 


For  the  general  initial-value  problem,  we  have  ro(0),  rx (0) , and 
rx(0).  From  these,  by  (2.4),  (2.7),  (2.42),  (2.9),  (2.43),  and  (2.47), 
we  may  calculate  V,  E,  t*,  TT^  , and  t*.  If  TTj  > 0,  there  will  be  a to  € (0,T*] 
such  that 


Ri(0)  = F(t0) 


When  rx (0)  > 0,  the  flow  for  Outfit*  (t*-to)  is 


ri  (t)  - W3  f(t0  + -^ 


(2.50a) 
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for  t*  (t*-t0)  * t < t*(2T*-T0)  it  is 


ri(t)  - W3  F (2t*-t0-  ^ 


and  for  t > t*  (2t*-to)  it  is 


rx(t)  = W3 


in  the  inelastic  case.  When  rj  (0)  < 0,  the  flow  for  0 £ t < t*T0  is 


ti(t)  = 


V^3  F 


(Jo-  ^*  ) , 


and  for  t > t*  tg  it  is 


rx(t)  « 


(2.50b) 


(2.50c) 


(2.51a) 


(2.51b) 


in  the  inelastic  case.  In  the  case  of  elastic  collisions,  these  flows  will 
be  time-periodic  with  period  2t*T*,  and  to(t)  will  vary  between  0 and 

**  CAT- 

At  this  juncture,  let  us  observe  that  one  or  the  other  of  the  free 
surfaces  will  tend  to  exhibit  Taylor  instability  at  some  point  during  the 
flow,  and  hence  the  solutions  with  non-zero  pressures  imposed  on  the  free 
surfaces  will  be  of  rather  limited  physical  interest.  It  follows  from  an 
analysis  elsewhere  (Ref.  3)  that  we  may  expect  Taylor  instability  at  the 
surface  r - r0(t)  if  and  only  if 


P|  < 0 

ar  r0(t) 


(2.52a) 
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and  Taylor  instability  at  the  surface  r = rx (t)  if  and  only  if 


£U<t>>0- 


(2.52b) 


Note  that  we  can  write 


P = 0 + Px 


(2.53a) 


r0(t)  rx(t) 


where  0 is  superharmonic  and  8 = 0 at  r0(t)  and  rx(t).  Hence, 


|^|  * 0 and 

3r  r0(t)  Sr'rx(t) 


s 0 


(2.53b) 


Accordingly,  from  (2.52)  and  (2.53),  we  will  never  get  Taylor  instability 
at  r = rQ(t)  when  Px  a 0,  or  at  r * r^t)  when  Px  £ 0. 

For  the  case  we  have  considered  so  far,  where  Px  > 0,  if  rx  (0)  > 0 then, 

, L / ^ \ ^ 

because  rx(t)  is  bounded  above  by  \r'  (1+  } , we  must  get  rx(t)  < 0 

\ ^ * 

for  some  t € (0,  t*(T*— r0)),  and  thus  Taylor  instability  at  r = rx(t).  An 
explicit  calculation,  based  on  (2.2),  (2.3),  (2.44),  and  the  known  values 
of  P at  r ■ r0(t)  and  r = rx(t),  yields 


V*/3 

E 


(2.54) 


Define  Rq  as  the  unique  positive  solution  of 
((r3°+1)i/j  - Rq)2 


_ 3 

1 2 I*3  //./n<0  ■ -‘1,S 


Ro  (4(Rq  +1)  - 3R?) 


+r,  3/3  . l/s  jO 

( (R6°+l)  + 2r£  (R*^  +1)  + 3Rq  ) . (2.55) 


i 


J 


I 

I 

I 


1. 
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It  is  easy  to  see  that  0 < R^  < ( ^ . For  r0(t)  > V1^3  R^,  we  get 

Taylor  instability  at  rx(t).  If  rx  (0)  < 0,  then  we  do  not  get  Taylor 
instability  at  the  boundary  r = rx  (t)  unless  ro(0)  > V1^3  Re- 
consider the  case  where  Px  < 0 but  e in  (2.42)  is  > 0.  We  may  still 
use  equation  (2.45).  Define 


+ i 6-  6-  ?)‘ 1 . 

7 =3  J,  7777 — 2 . 

1 a- i)1/a  (i-  f ttx  (§-i)j 


(2.56) 


•Jf  -J- 

Then  we  will  have  Vx  t “ as  t j t t.  This  follows  also  from  the  equation 
of  motion 


..  _ _ 1_  dPi  _ 

ri  " p0  dr  rx  " p0  V 


as  rx  | 00 


(2-57) 


For  0 < t < t , (2.45)  gives  Rx  as  a monotonically  increasing  function  of  t: 


Rx  = F(t)  . 


(2.58) 


In  the  same  way  we  derived  equation  (2.54),  we  find 


i 

JL 


1 

I 

1 


W3  i api 

E p0  dr 


to(t)’<i-5r)'n,+ 


(2.59) 


It  is  easy  to  see  that  there  is  a unique  positive  solution  Rq  of 
- "x  - | ■ (3(r£3+1)3*  + 2R?(rJ3+1)1/3+  r£3  ) . 


For  r0(t)  > W3  r£,  we  get  Taylor  instability  at  r0(t). 


(2.60) 
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Now  consider  the  general  initial-value  problem,  where  we  are  given  r0 (0), 
rx(0),  and  rx(0).  From  (2.4),  (2.7),  (2.42),  (2.9),  (2.43),  and  (2.56) 

'it 

we  compute  V,  E,  t , TTX  , and  t . We  find  to  € (0,t  ) as  the  solution  of 

Ri  = F(t0)  . (2.61) 

If  rx (0)  > 0,  the  flow  for  0 ^ t < t*  (t+-t0)  is  given  by 

ri  (t)  = W3  F(^o+^-  (2.62) 

For  t ^ t*  (t+-td)  the  flow  is  undefined.  If  rx (0)  < 0,  the  flow  for 
0 £ t < t*T0  is  given  by 


ri  (t)  = f'3  F^0-  f*)  • 

In  the  inelastic  case,  for  t > t*  T0  we  have 

ri  (t)  = V*7*  , 


(2.63a) 


(2.63b) 


JL  ^ 

and  in  the  elastic  case,  for  t td  < t < t (t0+t  ) we  have 

ri  (t)  = V173  F^  - ^ . (2.63c) 

i 

Suppose  finally  that  e as  given  by  (2.42)  is  < 0.  We  write  E = -e. 

Define  t*,  T,  Ro , Rx , or,  and  Vx  by  (2.9),  (2.10),  and  (2. lid),  and  let 


x*,.k.  mi. 
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„ _ V Pi 
1 E p0  • 


(2.64) 


The  equations  of  motion  in  dimensionless  form  are  given  by  (2.11a),  (2.11c), 


(i -*)- 


l - 


(2.65) 


It  is  clear  that  for  this  flow,  for  all  time. 


Ra  2 L. 

80  2^1  ' 


(2.66) 


Thus,  the  liquid  remains  away  from  the  origin,  and  we  do  not  have  the  possi- 
bility of  "collisions",  either  elastic  or  inelastic. 

3 

Consider  the  flow  with  V!  (0)  = - -r — - + 1: 


= 4 J 


J1  vH1-  t)  ) 

- 2?,+1  <?-1>l/a  (-1- 


(2.67) 


I 

1 

I 

i 

I 


We  take  positive  signs  for  the  square  roots.  (2.67)  gives  Rx  as  a monotone 
increasing  function  of  t: 


ri  - F(t)  . 


(2.68) 


As  t 1 0 we  get  from  (2.67) 


3 9 

v>  - 1 - 2*7  + f 


o-*r 


(2.69a) 
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**('-$* 

((i-*y -(-to 


(2.69b) 


Define  T by 


,..i  j (» if)'" 

3 - 2ST+1  <5-l)1/e 


(2.70) 


Using  (2.2),  (2.3),  (2.65),  and  the  prescribed  values  of  P at  the 
free  surfaces,  we  get 


V3 

v i_  apj 

E p0  dr1 


ro  (t) 


/ l+  i 2°  + 1 & 

3 »' IS. 


Ri  R? 


1 /'-L.  J_\  (2-  + 2 _i_\l 

" 2 VRo  “ RiArI  R0R1  R?/(  * 


(2.71) 


It  is  clear  from  (2.71)  that  we  get  Taylor  instability  at  the  surface 
r = r0(t)  throughout  this  flow. 

For  the  general  initial-value  problem,  we  are  given  ro(0),  ^(0),  and 
^(0).  Using  (2.4),  (2.7),  (2.42),  the  equation  E = - e,  (2.9),  (2.64),  and 
(2.70),  we  compute  V,  E,  t*,  TTj , and  t+.  We  find  To  € [0,t+)  as  the  solution 


Ri  - P(T0) 


(2.72) 
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If  rx (0)  > 0,  the  flow  for  0 £ t < t*  (t+-t0)  is 

MO  = V173  F (t0  + -^  . (2.73) 

The  flow  is  undefined  for  t £ t*  (t+-t0).  If  ^(0)  < 0,  for  0 £ t £ t*  t0  » 
we  have 

ri  (t)  = V173  F (r0  - ~A  (2.74a) 

and  for  t*  T0  £ t < t*  (to+t+)  » 

ri  (t)  = V173  fQ*  - • (2-74b> 

For  t ^ t*  (t0+t+)  the  flow  is  undefined. 

TWO-DIMENSIONAL  FLOW,  PRESSURE  APPLIED  AT  ONE  OF  THE  FREE  SURFACES 

The  velocity  potential  is  given  by  (2.23)  and  the  kinematical  free  sur- 
face conditions  by  (2.24).  We  set  the  pressure  on  the  inner  surface  equal  to 
0 and  assume  that  the  pressure  assumes  a value  Px  that  is  constant  over  the 
outer  surface  r = rx . The  conservation  of  volume  is  given  by  (2.25).  When 
Pi  is  constant  in  time,  the  conservation  of  energy  takes  the  form 

r,  (O  P,  rx  (0)  p, 

A2(t)  *n~7-T  + -i-  t?(t)  - A2(0)  In  — rg(0)  ■ c . (2.75) 

ro(c7  Po  ro(u7  Po 

e as  defined  here  will  be  > 0 whenever  Pi  > 0,  and  sometimes  when 
Pi  < 0.  When  e > 0,  we  set  E = e.  We  define  t*,  t,  Ro » Ri  , or,  and  Vi 


I 
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by  (2.27),  (2.28),  and  (2.29d).  Let 


n v £l 

1 = E Po  * 


(2.76) 


The  equations  governing  the  flow  are  (2.29a),  (2.29c),  and 


Ri 


a2  in  — = 1 - TTj  (Ri  -1) 


(2.77) 


Let  us  examine  the  special  flow  with  Vi  (0)  = 1 (Ro(0)  = 0) : 


3/3 


2/2  1 (l-nx(5-l)  ) 


1/8 


(2.78) 


As  usual,  we  take  positive  signs  for  all  square  roots.  The  behavior  for 
small  T is  given  by  (2.33a,b). 

Suppose  Px  > 0.  Then  it  is  clear  from  the  energy  equation  (2.77)  that 
Vx  will  have  an  upper  bound. 


- 1+S7  * 


(2.79) 


We  define  t by 


2/2  1 (l-TTi(S-l)) 


1/3 


(2.80) 
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and 


Ro 


. SSil 


•Va 


Xn(l-KTx  ) 


(T*“T 


(2.81b) 


For  0 < t < t*,  (2.78)  gives  Rx  as  a monotonically  increasing  function  of  t, 
which  we  denote  by 


Ri  = F(t)  . 


(2.82) 


As  we  have  seen  in  the  three-dimensional  case,  if  we  get  Taylor  instability 
at  all  it  will  occur  at  the  surface  r = rx(t).  Using  (2.2),  (2.23),  (2.77),  and 
the  prescribed  values  of  the  pressure  on  the  free  surfaces,  we  get 


W3  l_  dPi 

E p0  Sr'r^t) 


' ri)  | • 


(2.83) 


Ro 


We  can  check  that  there  is  only  one  positive  root  of  the  equation 


in 


- TT, 


(l+  ' iT 


in  ^1+ 


*) 


(2.84) 


+ 1 


* f \ V/2 

and  this  root  satisfies  0 < Ro  < ^ — j 


Now  let  us  consider  the  general  initial-value  problem,  where  we  are 
given  ro(0),  ^ (0) , and  ^(0).  Using  (2.24),  (2.25),  (2.75),  the  relation 
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E - €,  (2.27),  (2.76),  and  (2.80),  we  find  V,  E,  t*,  , and  t*.  There 

will  be  a t0  € (0,t*]  with 


Ri(0)  = F(t0  ) . 


If  ^ (0)  > 0,  the  flow  for  0 £ t £ t*  (t*-t0)  is 


MO 


'(T°  + ?) 


(2.85a) 


for  t*  (t*-t0)  s:  t < t*  (2t*-t0X  it  is 


ri  (t)  - W3  F ^2r*  - to-  ^ , 


and  for  t > t*  (2t*-t0)  it  is 


(2.85b) 


ri  (t)  = W3 


in  the  inelastic  case.  In  the  elastic  case,  for  t*  (2t*-t0)  < t £ t*  (3t*-t0)> 
we  get 


(2.85c) 


MO  - W3  F^~  - 2t*+t^ 


(2.85d) 


and  a flow  that  is  periodic  in  time  with  period  2t*  t*.  For  these  flows 
we  will  have  Taylor  instability  at  r « rx(t)  whenever  r0(t)  > W3  Ro»  which 
will  eventually  occur  if  rl (0)  > 0.  If  rt (0)  < 0,  the  flow  for  0 t < t*  t0  is 


Mt)  - V173  F ^r0  - ^ 


(2.86a) 
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and  for  t > t*  To  it  is 


ti  (t)  = W2 


(2.86b) 


in  the  inelastic  case.  In  the  elastic  case,  for  t*  t0  < t £ t*  (t0+t*)>  we 
have 


ri  (t)  - V1/2  - t0J  , 

and  for  t*  (to+t*)  ^ t < t*  (T0+2T*),we  have 

^(O-V^2  F^2t*+t0- 


(2.86c) 


(2.86d) 


the  flow  being  periodic  in  time  with  period  2t*  T*.  In  the  inelastic  case 
we  will  not  get  Taylor  instability  at  r = rj,(t)  if  rx  (0)  < 0 and  ro(0)  £ 

V1  ^ r£.  Otherwise.at  some  point  in  the  flow, Taylor  instability  will  occur 
on  the  outer  surface. 

Next  we  suppose  that  e in  (2.75)  is  > 0,  but  < 0.  Equation  (2.78) 
still  holds  and  gives  Rx  as  a monotonically  increasing  function  of  t: 


Ri  - F(t)  , 


(2.87) 


for  0 £ t < 08 . From  (2.78)  we  see  that,  when  Vx  f 

_ _1_  ( jA1/2  In  Vx  (2.88a) 

Zj2  ^ 

or 

V2  (-TTi)1'2  T 
Ri  - e 


(2.88b) 
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We  can  also  verify  this  directly  from  the  equation  of  motion, 


1 dPi 

r,  « - 

1 Po  3*  rx 


■ p‘  r>  as  r>  ’ • 


(2.89) 


As  we  have  seen,  Taylor  instability  can  only  occur  at  the  inner  boundary. 
By  using  (2.2),  (2.23),  (2.77),  and  the  prescribed  values  of  P at  the 
boundaries,  we  find  that 


V3^  1_  dPi 
E Po  Z*  r0(t) 


*°ln  R^ 


^ R^ 


(2.90) 


There  is  a unique  positive  solution,  Rq,  of  the  equation 


(1+£)  K1**-)'*""1 


(2.91) 


For  r0(t)  > V1/2  R^,  we  will  get  Taylor  instability  at  r * r0(t). 

In  the  general  initial-value  problem  we  are  given  r0 (0) , rx (0) , and 
rx(0).  We  use  (2.24),  (2.25),  (2.75),  the  equation  E - e,  (2.27),  and 
(2.76)  to  find  V,  E,  t*,  and  nx  . We  find  t0  € (0,®)  such  that 


Ri(0)  - F(t0)  . 


If  rx (0)  > 0,  the  solution  for  t a 0 is 


rx(t)  - F ^t0  +^\ 


(2.92) 


--MBi ■ 111  LX. I 111  i 
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31. 


If  rx(0)  < 0,  the  solution  for  0 £ t < t*  t0  is 


ri  (t)  “ F (r0  - pf) 


and  for  t > t To  it  is 


ri(t)  = W3 


(2.93a) 


(2.93b) 


in  the  inelastic  case.  In  the  elastic  case,  for  t > t*  To > we  have 


ti  (t)  = W2  f(%  - T0^  • 


(2.93c) 


There  will  be  no  Taylor  instability  in  the  inelastic  case  if  rx  (0)  < 0 and 
ro(0)  £ V1  2 R^.  Otherwise,  at  some  point  in  the  flow, we  will  get  Taylor 
instability  on  the  inner  surface. 

Now  suppose  e in  (2.75)  is  < 0.  We  write  E * - e.  We  use  (2.27),  (2.28), 
and  (2.29d)  to  define  t*,  T,  Rq , Rx , a,  and  Vx  , and  we  let 


-I  r 

E p0 


(2.94) 


In  dimensionless  variables,  the  equations  of  motion  are  (2.29a),  (2.29c),  and 


a2  la  ~ = - 1 - tti(r?-1)  . 


(2.95) 


We  must  have  R§  & - , and  thus  we  cannot  have  "collisions"  of  liquid  at 

the  origin. 
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32. 


Let  us  examine  the  special  flow  with  Vx (0)  = 1 - . We  get 

t._l 


Z/2  1_  (-1-^(5- 1)) 

TT, 


. 


(2.96) 


and  taking  positive  signs  for  the  square  roots,  this  gives  Rx  as  an  increas- 
ing function  of  t € [0,“): 


Ri  - F(t)  • 

The  asymptotic  behavior  as  T t 00  is  given  by  equation  (2.88). 
The  behavior  as  t I 0 can  be  obtained  from  (2.96): 


Vi 


1 - 


2TT, 


TTj  £n(  1-tTj) 


and 


Ro 


L i.V/s  , <~n>>3/3  Ta 

\ "i/  *n(l-TTi) 


(2.97) 


(2.98a) 


(2.98b) 


To  study  the  stability  of  the  surface  r = r0(t),  we  use  (2.2),  (2.23), 
(2.95)  and  the  prescribed  values  of  F at  the  boundaries  to  get 


y®  fa 

E 


1 dPl 

* S? 


r0  (t) 


Ro  An 


*•  1 A 


(2.99) 


and  this  is  < 0 throughout  the  flow.  Thus,  we  geL  Taylor  instability  at 
t - ro(t). 

To  solve  the  general  initial-value  problem,  we  suppose  that  we  are 
given  ro(0),  rx  (0) , and  rx(0).  We  use  (2.24),  (2.25),  (2.75),  the  equation 
E ■ - c,  (2.27),  and  (2.94)  to  compute  V,  E,  t*,  and  nx  . We  find  t0  € [0,®) 


- — - — 
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so  that 


I k(0)  - P(T0)  . 


Then, if  rx (0)  > 0,  the  flow  for  t ^ 0 is  given  by 


Cl  (t)  = ^ F (r0  + ^ 


(2.100) 


If  rx  (0)  < 0,  the  flow  for  0 £ t £ t T0  is 


and  for  t ^ t t0 


ri  (t)  - Wa  P (to  - “f) 

rx(t)  - W2  - t0^  . 


(2.101a) 


(2.101b) 


TWO-DIMENSIONAL  FLOW  WITH  CIRCULATION,  ZERO  PRESSURE  ON  THE  FREE  SURFACES 


In  place  of  the  velocity  potential  (2.23),  we  have 


9 - A(t)  In  + C 6 + B(t)  , 


(2.102) 


where  8 is  the  azimuthal  angle.  The  important  case  is  where  C is  independent 
of  time  as  well  as  space.  The  kinematical  boundary  conditions  are  still  given 
by  (2.24).  The  dynamical  boundary  conditions  are  obtained  from  (2.2),  (2.102), 
and  the  fact  that  P - 0 on  the  free  boundaries: 


. r°  A 1 .1.1  A3+(f 

A in- — rx  + B + 7T  —5 — - 0 , 

ri  ri  *■  rj 


(2.103a) 


A . , j,  , 1 A S+Ca 

‘ 7T  1 2 — 0 • 


(2.103b) 
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34. 


From  the  kinematical  conditions,  conservation  of  volume  follows  as 
in  (2.25).  For  the  conservation  of  energy,  we  get 


(A3(t)+C3)  in  - (A3(0)+C2)  in  = E . 


(2.104) 


We  may  use  (2.27),  (2.28),  and  (2.29d)  to  define  t , t,  Rq  , Ra  , or, 
and  Vx . Further,  we  let 


(2.105) 


The  governing  equations  are  (2.29a),  (2.29c),  and 


(a2  +Y2 ) in  ~*  = 1 . 
«o 


(2.106) 


Since  Y is  constant,  it  is  clear  from  (2.106)  that  we  cannot  have  Rq  * 0. 
Thus,  these  flows  will  never  exhibit  "collisions".  The  minimum  possible 
value  of  Rq  will  be  achieved  when  r0  = A = a = 0: 


Y2  in  — — « 1 *»  Vi  = • -* 

1 l-e-2/YS 


(2.107) 


Let  us  consider  the  flow  for  which  V*  (0)  is  given  by  (2.107).  Then  we 


2*/2  L-.V'  V *n(l+  fT ^ 

l-e"2^  V V 


(2.108) 
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In  this  equation  we  take  the  positive  sign  for  the  square  root,  and  thus 
we  get  Rx  as  a monotonically  increasing  function  of  T: 


Ri  - F(t) 


From  (2.108)  we  may  obtain  the  behavior  as  T t 0; 


(2.109) 


v,-— 1 — + , 

. -2/v2  2 \ J 


i/i  i vV3  i ( i i 

2v^  V Y5  ~ Y1 1 V*  I r3  Y1 
Ro  - e T \ eT  -e  X j +±~e  "e  / 


(2.110a) 


(2.110b) 


The  asymptotic  behavior  as  t t 00  is  given  by  (2.34). 

Now,  suppose  we  have  r0 (0) , rx (0) , and  £*(0).  From  (2.24),  (2.25), 
(2.104),  (2.105),  and  (2.27)  we  find  V,  E,  Y and  t*.  We  find  t0  € [0,°°) 


such  that 


Ri(0)  = F(T0) 


If  rx (0)  > 0,  the  solution  for  t £ 0 is  given  by 


rx  (t)  = \A2  F 


(To  + 


(2.111) 


If  rx (0)  < 0,  the  solution  for  0 s t ^ t t0  is 


and  for  t £ t*  To 


ri  (t)  = V1/2  F^t0  - 


ri(t)-vV*  - T0j  . 


(2.112a) 


(2.112b) 
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36. 


TWO-DIMENSIONAL  FLOW  WITH  CIRCULATION,  AND  PRESSURE  APPLIED  AT  ONE 
OF  THE  FREE  SURFACES 


We  treat  now  the  case  where  there  is  circulation,  ar.d,  as  well,  a 
pressure  P0,  which  is  constant  over  the  surface,  is  applied  at  the  inner 
boundary.  The  pressure  at  the  outer  boundary  is  taken  to  be  0. 

The  velocity  potential  is  given  by  (2.102).  The  kinematical  boundary 
conditions  and  conservation  of  volume  are  given  by  (2.24)  and  (2.25).  For 
the  dynamical  boundary  conditions,  we  insert  (2.102)  into  the  Bernoulli 
equation  (2.2)  and  get 


+ B 


♦* 


A2+C2 


po  o 
+ — = 0 
Po 


(2.113a) 


A . . 1 A 2+C3 

- ~ + B + = 0 * (2.113b) 

ri  1 r? 

The  expression  of  conservation  of  energy  assumes  the  form 

(A2(t)+C2)  In  ^ r|(t)  = (A2  (O)-t-C2 ) In  r2  (0)  = e . (2.114) 

rOKt)  Po  ro(0)  Po 


If  P0  < 0,  e will  be  > 0.  Ihere  will  also  be  cases  when  e > 0 and 
PQ  > 0.  So  first  we  consider  e > 0 and  set  E = e.  We  define  t*,  t,  Ro  » 
Rt , a,  Y,  and  Vx  by  (2.27),  (2.28),  (2.105),  and  (2.29d).  Let 


"o 


V Po 
E po 


(2.115) 


The  equations  governing  the  flow  are  (2.29a),  (2.29c),  and 


(or2  +Y2 ) Xn  “ 1 + tt0(R?-1)  . 
«o 


(2.116) 
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37. 


Since  y is  constant,  we  cannot  have  Ro  = 0,  and  there  will  be  no  "collisions". 
The  minimum  value  of  Ro  will  be  achieved  when  Of  = 0: 


V2  t n — = 1 + TT0  R§ 
Ro 


(2.117) 


When  ttq  s 0 there  will  be  one  and  only  one  root  Ro  > 0 of  (2.117).  We 
call  this  root  Ro(Y»tto)a  For  this  case*  the  flow  will  have  rD(t)  £ V1^2  Ro  (Y^o) 
Given  Y and  tt0  < 0,  there  may  be  0,  1,  or  2 roots  Rq  > 0 of  (2.117).  If 
there  are  no  roots, 

Ri 

Y2  l n — > 1 + tt0  r§ 

Ro 

for  all  Rq  > 0.  On  the  other  hand,  from  (2.116),  for  any  real  flow 


Y2  in  ^ 1 + n0  R§  , 
Ro 


and  thus  for  such  flows  equation  (2.117)  will  have  1 or  2 roots.  We  easily 
check  that  if  (2.117)  has  just  one  root,  we  must  have 


in  (\+  — — ^ 
Y R|(0)/ 


= 1 + TT0  R§  (0) 


and 


Y2 


2Rq  (0)  R?(0) 


= tt, 


o * 


from  which  we  get 


(2.118a) 


(2.118b) 


1 

li 


A(0)  = 0 


and,  on  combining  (2,113a)  and  (2.113b), 


A(0)  = 0 


(2.119a) 


(2.119b) 


r 


nm  i»i|ni,».'  i u.n» ,) 


T 
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38. 

In  this  case,  the  flow  is  steady,  and  the  negative  pressure  P0  on  the 
inner  boundary  is  just  sufficient  to  balance  the  centrifugal  force  of  the 
swirling  motion.  So  we  may  consider  the  case  where  (2.117)  has  two 
distinct  positive  roots.  We  call  them  Ro1^(V>TTo)  and  RoS^(V»TTo)»  with 
Ro^the  smaller  root.  Thus,  the  flow  will  satisfy 

^^(V.TTo)  ^ r0(t)  £ V^a  4a)(Y,TT0)  (2.120) 


for  all  t ^ 0. 

A more  complete  description  of  the  restrictions  imposed  on  -TT0  and 
V2  by  the  requirement  that  (2.117)  have  at  least  one  solution  is  as  follows: 

We  must  have,  for  0 < - TT0  < ® and  0 < Y2  < °°, 

0 < Y2  S G(-tt0)  (2.121a) 

and 

0 < - TT0  <.  G_1  (Y2 ) , (2.121b) 


where  the  curve  Y2  = G(-TT0)  or  - TT0  = G 1(Y2)  is  given  parametrically  by  the 
equations 

- "o  = — 7 =T  (2.122a) 

§ [(1+5)  In  (l+  + 1^ 

and 


(2.122b) 


for  0 < 5 < ®.  We  see  that 

\ 

^ <-"o)  < 0, 

\ 

\ 

. 1 \M 
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39. 


and  hence 


<4i2l<0,  i c-  M < 0,  0<K<». 

dx  dx 


(2.123) 


From  (2.122)  we  obtain 


and 


also 


and 


-l 


1 

* 2x 

as 

x l 0 

(2.124a) 

2 

in  x 

as 

x t “ ; 

(2.124b) 

x 2/x 
2 e 

as 

x i 0 

(2.125a) 

1 

"*  2x 

as 

x T “ . 

(2.125b) 

Now  for  tt0  c 0,  let  us  consider  the  special  flow  with  Vx  (0) 
We  get 


(Ro 


o(1))%  1. 


i r 

T = — J 

2^2  ,D(i)  2 


in 


(Rov^)  +1  \1  + tto(S-1)  - f-  in(l+ 


d§ 


(2.126) 


Taking  the  positive  square  root,  we  see  that  (2.126)  gives  Sj  as  a monotonically 


increasing  function  of  t for  0 < t < t , 

Ri  * F(t)  , 


(2.127) 


where 


,*  1 


(Ro(2))%1 


I 


ini 


0*  fh) 


j/a 


^ (Ro(i))%i  v + w-1*  - £ H1+ 


d§  . (2.128) 
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40. 


As  t 1 0,  we  have 


Vx  - (Ro1*)2  + 1 + 


Y3  . Y2  | 

( 1 

1 \ 

1 + TT0  <*<>>)».  2 

^(Rq1^)2 

(Ro  1^)2+1  _ 

Ta  (2.129a) 


and 


Ro  -•  Rq  ^ + 


2R^l)(l-HT0(R0v 


0(l))a)  _ 2 Cro(i))2  (Ro^^+li 


(2.129b) 


At  the  other  end  of  the  range,  we  find  the  behavior  as  (t*-t)  i 0: 

_ 

Vx  - (Ro8))a  + 1 ^v\-  tt0  - (—A A ) (t*-t)2 

1-HT0(Ro(3))2  . 2 '(R^)2  (Rq2^)2+1  — (2. 


130a) 


and 


Ro  - Ro 


(a)  _ 


2r£ 


y2 r_  n . t_( i i. ..  s 

V-HTqCRo^)2)  L 0 2 MRoa))2  (Roa))a+l/. 


(T*.T)a.  (2.130b) 


If  we  get  Taylor  instability,  it  will  occur  at  the  surface  r = rx(t). 

From  (2.2),  (2.102),  (2.116),  and  the  known  values  of  P on  the  free  boundaries, 
we  find 


A3  1 3Pi  1 ( 

E P°  8t,tl(t)\in| 


_ 1 

0+[l-HT0(^-l)][  - ■**>  " ^ 

Rf  . Ri 
ln  R§ 


(2.131) 


As  in  the  case  where  tnere  was  no  circulation,  there  is  a number;  Rq,  dependent 
* a 1 

only  on  -tt0  , 0 < R©  < - — , such  that  the  expression  in  (2.131)  is 

T1o 

* * * 
negative  for  Rq  < Rq  and  positive  for  Rq  > Rq.  R<,  satisfies 


in 


TTn 


A & 


(2.132) 


+ 1 


M 


L 

r 


_LL 
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41. 


The  question  of  Taylor  instability  thus  depends  on  the  relation  of 


l£(-Tr0)  to  8o1^(V»^o)  and  Ro  **  (V.^o) • We  see  from  (2.117)  that  as  Vs  4 0 


,(») 


we  get 


®oX\v»n0)  - e 2y3^e^ - e ^ 


_ l/s 


) 


- 0 , 


(2.133) 


as  in  (2.110b).  And 


R,CaW>  - • 


(2.134) 


We  recall  that  the  expression 

f-  in  (l+  - 1 - n0  § 

vanishes  when  5 is  (Fq1^)3  or  (Ro8^)8  and  is  negative  for  (Rq1^)2  < 5 < 


(Ro(2))2.  Hence,  differentiating  with  respect  to  § at  (R^1^)2  and  (Rq^)2* 


we  get  that 


25(1+5) 


- nr 


is  negative  at  (Ro1^)8  and  positive  at  (Ro^)2*  Now  differentiate  (2.117) 
with  respect  to  y2  when  Rq  is  R©1^  or  R^8}  We  find 


a*Sl)<Y,nb) 


dY 


2 0 


(2.135a) 


and 


aRo(S)(Y,TT0) 


dY 


s 0 


(2.135b) 
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Sevas  Y3  increases  from  0 to  G(JT0).  the  interval  f Rq1  ^ , R<$3^  shrinks 

f / i \Va\  ^ _ ' 

from  fO,  (-  — J j to  a small  neighborhood  of  So,  where  Rq  satisfies 
(2.118)  with  Y3  replaced  by  G(-n0).  it  is  easy  to  see  that 


So  > • (2.136) 

For,  as  we  saw  in  (2.119),  the  corresponding  limiting  flow  is  stationary, 
and  the  only  acceleration  of  fluid  elements  at  r - V1^2  (So+1)1^3  is  the 
centripetal  acceleration  inward,  which  can  only  be  caused  by  a positive 
value  of  ^ at  that  point.  As  Y3  increases  from  0 to  G(-ttq),  it  will 
cross  a value  H(-ttq)  such  that 

( (H(-tt0))i/3,tTo)  - Ro  • (2.137) 

For  0 < Y3  < H(-tt0),  the  surface  r ■ rx(t)  will  be  stable  when  W3  Ro1^(Y,T1o) 

*•  rQ(t)  s V^3  Rq  «nd  unstable  when  V^3  Ro  < r0(t)  £ V^Rq (Y.^o)  • for 
H(-n0)  < y3  < G(-n0),  the  surface  r - rx(t)  will  be  unstable  throughout 
the  motion.  He  may  compute  H(-tt0)  from  (2.137),  (2.132),  and  (2.117).  The 
result  is 

H - 2 ( l-n0 ) R*a  . (2.138) 

For  the  general  initial-value  problem,  with  P0  < 0,  let  us  suppose  ve 
have  given  ro(0),  rx (0) , and  ri(0).  He  use  (2.24),  (2.2S),  (2.114),  the 
equation  E - «,  (2.105),  (2.27),  (2.115),  and  (2.128)  to  find  V,  E,  Y,  t*, 

TT0,  and  t*.  Using  (2.127),  we  find  r0  € [0,T*1  such  that 


Ri  (0)  - F(t0) 


r 


L 
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If  rx  (0)  * 0,  the  flow  for  0 £ t £ t*  (t*-to)  is  given  by 


rl  (fc)  - ^3  f(t0  + pp) 


for  t*(T*-t0)  £ t £ t*  <2t*-t0),  by 


rl  (t)  « W2  F ^T*-To-  pp)  ; 


for  t (2t*-t0)  £ t £ t*  (3t*-t0),  by 


and  in  general,  by 


ri  (t)  - V172  F ^-pr  - 2t*+t0)  ; 


ri  <fc)  * ri  (t+2t*  t*)  . 


If  ix (0)  £ 0,  for  0 £ t £ t*  T0  we  have 


Mt)  - V*AF^0  - ; 


for  t*T0  £ t £ t*  (t0+t*). 


ri  (t)  - V172^  - to)  ; 


for  t*  (t0+t*)  £ t £ t*  <t0+2t*). 


fi(t)  - V1^  F (2t*  + To  - pp)  ; 


43. 


(2.139a) 


(2.139b) 


(2.139c) 


(2. 139d) 


(2.140a) 


(2.140b) 


(2.140c) 


and,  in  general,  we  have  the  periodicity  (2.139d). 
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44. 


Suppose  now  that  e given  by  (2.114)  is  > 0 but  PQ  > 0.  As  we  have 
seen,  there  is  one  positive  root  RoW^o)  (2.117).  Consider  the  special 
flow  with  Vx(0)  - (Ro(Y,TT0))a  + 1.  Then 


JL  }'  ( X d5 

^ (Ro(Y,TT0))a+l  ^0(5-1)-  f-  in(^1+  5^l)/ 


(2.141) 


If  we  take  the  positive  square  root  here,  we  get  Rx  as  a monotonically  increasing 
function  of  t for  t > 0: 


Ri  - F(t)  . 


(2.142) 


The  behavior  as  t l 0 is 


Y H tt0+ 


Vx  - (®o(Y  ,Tro))a+l  + 


2 (Ro  (Y  ,"o  ) )3  ( (Ro  (Y  ,"o  ) )3+l). 
1 +tt0(Ro(Y,tt0))2 


(2.143a) 


Ro  - Ro (Y»TT0)  ^-A0  +2^<V»TTQ))3.W(Y>TT0))a-H))Ta  # 

2Ro(Y,TT0)(l-m0(Ro(Y,TT0))2) 


(2.143b) 


As  t t 00  we  get 


•JxQ 


Rx  -*  e 


(2.144) 


Regarding  Taylor  instability,  it  will  occur  only  at  the  inner  boundary, 
r ” r0 (t) . From  (2.2),  (2.102),  (2.116),  and  the  known  boundary  values  of 
P,  we  find 
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v3/3  i api  i 

E P°  Sr  ro(t)  R^n 


i ( U 

— j.no  + (l«0  («-!»(- — 


. *L 
ln H 


(2.145) 


As  In  the  case  without  circulation,  we  find  a unique  Ro » dependent  only 
on  ttq  , wi  th  0 < < “ and 


(2.146) 


When  rQ(t)  £ W3  Ro.  we  get  no  Taylor  Instability.  When  r0(t)  > V1^3  Ro, 
we  get  Taylor  instability  at  r = r0(t). 

So  the  only  question  is  the  relation  of  R^  to  the  lower  limit  Ro(Y,TTo) 
of  the  inner  dimensionless  radius.  It  is  easy  to  see  from  (2.117)  that  we 


Ro(Y,"o)  I 0 as  Y3  JO 


(2.147a) 


Ro(Y,tt0)  t » as  Y3  t 


(2.147b) 


As  Y2  increases  from  0 to  ®,  it  will  cross  a value  K(n0)  such  that 


Ro(K(n0),TT0)  - Ro  . 


(2.148) 


From  (2.148),  (2.117),  and  (2.146),  we  get 


K - 2(l+Ro"3)  . 


(2.149) 


When  Y3  > K(tt0),  the  surface  r ■ to(t)  will  be  unstable  during  the  whole 
motion. 
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Now  consider  the  general  Initial-value  problem,  where  we  have  ro(0), 
rl  (0) » and  ri(0).  From  (2.24),  (2.25),  (2.114),  the  equation  E = e, 
(2.27),  (2.105),  and  (2.115)  we  compute  V,  E,  t*,  Y,  and  tt0  , We  find 


To  € [0,®)  such  that 


Ri(0)=F(to)  . 


If  fj (0)  2 0,  the  flow  for  t ^ 0 is  given  by 


ri  (t)  = V1^  F (to  + “5F  J 


(2.150) 


If  rx(0)  £ 0,  the  flow  for  0 s t £ t*  T0  is 


and  for  t i t t0 


ri  (t)  = V1/3  fTt0  " -C) 


Ci  (t)  = V173  F - To') 


(2.151a) 


(2.151b) 


Finally,  we  consider  the  case  where  e as  given  by  (2.114)  is  < 0. 
Write  E = -€.  Then  we  use  (2.27),  (2.28),  (2.29d),  and  (2.105)  to  find  t*. 


t,  Rq>  Ri » or,  Vi  and  y and  we  set 


V po 


n°  " E Po 


(2.152) 


The  governing  equations  are  (2.29a),  (2.29c),  and 


(oa4Ya)  In  ~ - - 1 +TT0(^-l)  . 


(2.153) 
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For  this  flow  we  must  have  R©  a So(Y,n0X  where  Bo  is  the  unique  positive 
solution  of 

« - 1 + n0  ^ . (2.154) 

when  y8  10,  (2.155a) 

So  t “ when  y2  t ® , (2.155b) 


and 


dSo 

d7 


2 0 . 


Consider  the  flow  with  Vi (0)  - (So(Y,n0))3  + 1.  We  get 


T 


(2.155c) 


(2.156) 


This  equation,  with  positive  square  roots,  determines  Rx  as  a monotonically 
increasing  function  of  T for  T € (0,®): 


Ri  - F(t)  . 


As  t l 0 we  get 


Vx  - Bg  + 1 + 


2g^(^+l) 


•) 


/ 


1 \ 


(2.157) 


(2.158a) 
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Ro 


TTn  + 


So  + 


So  fn 


2*?  <*?+!) 

er 


(2.158b) 


As  T t 00  we  have  equation  (2.144). 

The  Taylor  instability  at  the  surface  r " r0(t)  is  checked  by  using 
(2.2),  (2.102),  (2.153),  and  the  given  values  of  P at  the  inner  and  outer 
surfaces: 


V*/3  1 3Pi  1 

f - "o  I 

i . dr.j?  1 ) 

E Po  Sr  r0(t) 

Ro 

1 

‘arhB 

c 

« in£  : 

(2.159) 


and  this  is  negative,  giving  us  instability  at  r = r0(t),  throughout  the 
flow. 

For  the  general  initial-value  problem  we  use  ro(0),  rx  (0) , and  rx  (0) 
to  compute  V,  E,  t*,  V,  and  TT0  from  (2.24),  (2.25),  (2.114),  the  equation 
E - -e,  (2.27),  (2.105),  and  (2.152).  From  (2.157)  we  find  t0  € [0,®)  such 
that 

Ri(0)  - F(t0)  • 

If  r^  (0)  2:  0,  the  flow  for  t * 0 is 

rx(t)  - W2  F (jo  + pfj  . (2.160) 

If  (0)  s 0,  the  flow  for  0 s t £ t*  t0  is 

rx(t)  - Wa  F(j0  - , (2.161a) 

and  for  t 2 t*  To » 

rx  (t)  - V1/3  ffe  - t0 J . (2.161b) 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

Laurel  Maryland 


49. 


n 


REFERENCES 

li 


.. 

li 
* » 

T ■» 

1! 

n 

\ * 
4* 

u 


1.  Joel  C.  W.  Rogers,  "Relation  of  Turbulence  to  Energy  Conservation, 
Hadamard  Instabilities,  and  Stochastic  Processes",  fourth  in  a series 
of  reports  which  will  constitute  part  of  the  paper  "Algorithms  for 
Hyperbolic  and  Hydrodynamic  Free  Boundary  Problems". 

2.  Joel  C.  W.  Rogers,  "Algorithm  for  Water  Waves,  Error  Bound  for  Pressure 
Computation",  second  in  series  of  reports  referred  to  in  previous 
reference. 

3.  Joel  C.  W.  Rogers,  "Water  Waves:  Analytic  Solutions,  Uniqueness  and 
Continuous  Dependence  on  the  Data",  Naval  Ordnance  Laboratory 
NSWC/WOL/TR  75-43  (1975). 


51 


INITIAL  DISTRIBUTION  EXTERNAL  TO  THE  APPLIED  PHYSICS  LABORATORY* 

Tha  work  reported  In  TG  1307  vu  don*  under  Navy  Contract  N00017-72-C-4401.  Thia  work  1*  rslsted  to  T*»k  ZL50,  which 
1*  mpportfd  by  the  Office  of  Naval  Reaearch  (Cod*  438). 


ORGANIZATION 


LOCATION 


ATTENTION 


DEPARTMENT  OF  DEFENSE 


Depertwent  of  the  Nev 


NAVPRO 

NAVSEASYSCOM 


DTNSRDC,  Annapolla  Leb 
DTNSRDC 


U.S.  Naval  Academy 

Office  of  Naval  Rea. , Branch  Office 


Office  of  Naval  Reaearch,  Scientific 
Lleieon  Group,  Amer.  Embaaay 
Office  of  Naval  Reaeerch 


NAVAIRSYSCOM 


Puget  Sound  Naval  Shipyard 
Charleaton  Navel  Shipyard 
Long  Beach  Neval  Shipyard 
Philadelphia  Naval  Shipyard 
Porteaouth  Neval  Shipyard 
Norfolk  Navel  Shipyard 
Pearl  Harbor  Neval  Shlpyerd 
Hunter*  Point  Neval  Shipyard 
Rare  talend  Nevel  Shipyard 
Naval  Underaea  Center 


Naval  Reaearch  Lab. 


Neval  Weapon*  Center 
Naval  Surface  Weapons  Center 


Alexandria,  VA  22314 


Laurel,  MD  20810 
Washington,  DC  20362 


Annapolla,  MD  21402 
Betheada,  HD  20084 


Annapolla,  MD  21402 

Boaton,  MA  02210 
Chicago,  IL  60605 
New  York,  NY  10003 
Paaadena,  CA  91101 
San  Franclaco,  CA  94102 

AP0  San  Franclaco,  CA  96503 
Arlington,  VA  22217 


Waahlngton,  DC  20361 


Washington,  DC  20376 

Bremerton,  VA  98314 
Charleston,  SC  29408 
Long  Beach,  CA  90801 
Philadelphia,  PA  19112 
Portsmouth,  NH  03801 
Portsmouth,  VA  23709 
FP0  San  Franclaco,  CA  96610 
San  Franclaco,  CA  94135 
Valleno,  CA  94592 
San  Diego,  CA  92132 


Waahlngton,  DC  20375 


China  Lake,  CA  93555 
Dahlgren,  VA  22418 


SEA-0322 
SEA-032 
SEA-03B 
SEA-03512 
SEA-09G3 
Tech.  Lib. 
Code  01 
Code  15 
Code  18 
Code  1521 
Code  1541 
Code  1552 
Code  1843 
Code  5641 
B.  Johnson 
Lib. 

Dir. 

Dir. 

Dir. 

Dir. 

Dir. 


Code  438 
Code  200 
Code  210 
Code  211 
Code  212 
Code  221 
Cod  s 473 
Code  480 
Code  481 
Code  102 IP 
AIR-03 
AIR-03B 
AIR-310 
AIR-5301 
AIR-50174 
Dir. 

SP-2022 

Conmandor 

Comander 

Commander 

Tech.  Lib. 

Commander 

Commander 

Lib. 

Tech.  Lib. 

Tech.  Lib. 

Code  2501 
Tech.  Lib. 

Code  4007 
Code  2627 
Code  4000 
Code  7706 
Code  8'  41 
Lib. 

Tech.  Lib. 

Comp.  & Anely.  Lab. 


Requests  for  copies  of  this  report  from  DoD  activities  and  contractors  should  be  directed  to  DOC. 
Cameron  Station,  Alexandria,  Virginia  22314  using  DOC  Form  1 and,  if  necessary,  DDC  Form  55. 


'Initial  distribution  of  this  document  within  the  Applied  Physio  Laboratory  has  bean  mad*  in  accordance  with  a list  on  fit*  in  th*  APL  Technical  Publication*  Group. 


52 


INITIAL  DISTRIBUTION  EXTERNAL  TO  THE  APPLIED  PHYSICS  LABORATORY* 


ORGANIZATION 

LOCATION 

ATTENTION 

Department  of  the  Now  (cont'd) 

Navel  Surf ece  Weapone  Center  (cont 'd) 

Silver  Spring,  MD  20910 

J.  Enig 

Lib. 

J.  Rogers 

Naval  Ship  Engineering  Center 

Philadelphia,  PA  19112 

Tech.  Lib. 

Hyattsville,  MD  20782 

Code  6034 

Code  6101E 

Code  6110 

Code  6114 

Code  6136 

Code  6140 

Navel  Poetgraduete  School 

Monterey,  CA  93940 

Lib. 

Naval  Undetveter  Syateas  Center 

Newport,  RI  02840 

Tech.  Lib. 

Naval  Coeetal  System  Lab. 

Panama  City,  FL  32401 

Tech.  Lib. 

Naval  Mleelle  Center 

Ft.  Mugu,  CA  93  041 

Tech.  Lib. 

Occanogrepher  of  the  Navy 

Alexandria,  VA  22332 

Naval  Oceenogrephlc  Office 

Washington,  DC  20373 

Commander 

Commandant  of  the  Marine  Corpe, 

Scientific  Advleor 

Washington , DC  20380 

Code  AX 

Department  of  the  Army 

Rea.  & Tech.  Div. , Army  Eng.  Reactora 

Croup 

Ft.  Bel voir,  VA  22060 

Army  Mobility  Equip.  Rea.  Center 

Ft.  Belvoir,  VA  22060 

Tech.  Doc.  Center 

Army  Mleelle  Command,  Red  atone 

Scientific  Info.  Center 

Redstone  Arsenal,  AL  35809 

Chief,  Doc.  Sec. 

Army  Reeearch  Office 

Research  Triangle  Pk. , NC  27709 

Office  of  Chief  of  Staff 

Washington,  DC  20310 

Chief,  Rea.  & Dev. 

Department  of  the  Air  Force 

Office  of  Scientific  Reaearch/NA 

Washington,  DC  20332 

APDRD-AS/M 

Washington,  DC  20330 

DEPARTMENT  OF  TRANSPORTATION 

Coeet  Guard  Hq. 

Weehlngton,  DC  20591 

Lib.  Sta.  5-2 

U.s.  GOVERNMENT  AGENCIES 

National  Aeronautics  and  Space  Adnln. 

1ASA  Scientific  and  Tech.  Info.  Fee. 

Beltimore/Wsshlngton  International 
Airport 

9ASA  Levi*  Reeeerch  Center 

Cleveland,  OH  44871 

W.  J.  Anderson 

CONTRACTORS 

tape*  Corp. 

Redwood  City,  CA  94063 

A.  Eahel 

Aeroneutlcel  Res.  of  Princeton,  Inc. 

Princeton,  NJ  08540 

CALSPAN  Corp. 

Buffelo,  NY  14221 

Dreper  Lab. 

Cambridge,  MA  02139 

Lib. 

Flow  Reeeerch,  Inc. 

Cambridge,  MA  02142 

S.  Orszag 

General  Technicel  Servlcee,  Inc. 

Kent,  WA  98031 

Loe  Angelee,  CA  90045 

D.R.S.  Ko 

Teedon,  PA  19050 

A.  S.  Iberall 

Hydroneutlce,  Inc. 

Leurel,  MD  20810 

V. Johnson 

M.  P.  Tulin 

JATCOR 

Del  Her,  CA  92014 

R.  Chan 

J.  A.  Young  j 

Oceanic*,  Inc. 

Pleinviev,  NY  11803 

P.  Kaplan 

Poeeldon  Raeeerch 

Loe  Angeles,  CA  90049 

S.  Crow 

Shaker  Reteerch  Corp. 

Belleton  Lake,  NY  12019 

C.  Pan 

Southveet  Reeeerch  Inetltute 

Sen  Antonio,  TX  78228 

H.  N.  Abremson 

Ed.,  App.  Mech.  Rev. 

Tetre  Tech,  Inc. 

Peaadena,  CA  91107 

R.  Wede 

TRW  System.  Group 

Redondo  Beach,  CA  90278 

E.  Beum 

Attribution  of  tbit  docwmnt  within  tht  Appliwi  Physics  Laboratory  hot  boon  mad.  in  accordance  with  o Ktt  on  file  in  tho  APL  Tocbnicol  bublicotiont  Group. 


f 

a, 

- 


53. 

INITIAL  DISTRIBUTION  EXTERNAL  TO  THE  APPLIED  PHYSICS  LABORATORY* 


ORGANIZATION 

LOCATION 

ATTENTION 

No.  of 
Copwt 

UNIVERSITIES 

Brown  University 

Providence,  RI  02912 

J.  Clarke 

1 

J.  Liu 

1 

California  Institute  of  Technology 

Pasadena,  CA  91109 

H.  W.  Llepnann 

1 

M.  S.  Pleaaet 

1 

A.  Roahko 

1 

T.  Y.  Wu 

1 

Carnegie  Inatitute  of  Technology 

Pittsburgh,  PA  15213 

R.  C.  MacCamy 

1 

Case  Western  Reserve  University 

Cleveland*  OH  44106 

E.  Reohotko 

1 

Columbia  University 

New  York,  NY  10027 

V.  Caatelll 

1 

H.  G.  Elrod 

1 

Cornell  University 

Ithaca,  NY  14851 

E.  L.  Realer 

1 

Harvard  University 

Cambridge*  MA  02138 

G.  BIrkhoff 

1 

G.  F.  Carrier 

1 

Loa  Alamos  Scientific  Lab.,  Univ.  of 

California 

Loa  Alamos,  NM  87544 

C.  W.  HIrt 

1 

Massachusetts  Inatitute  of  Tech. 

Cambridge,  MA  02139 

M.  A.  Abkovitz 

1 

Com.  Off.,  NROTC 

Nav.  Admin.  Unit 

1 

L.  N.  Howard 

1 

P.  Mandel 

1 

C.  C.  Mei 

1 

E.  Mollo-Chrlstensen 

1 

J.  N.  Newman 

1 

R.  F.  Probatein 

1 

Mlaalaalppl  State  University 

State  College,  MS  39762 

J.  Thompaon 

1 

New  York  University 

New  York,  NY  10003 

J.  J.  Stoker 

1 

Pennsylvania  State  University 

University  Park,  PA  16802 

J.  L.  Lumley 

1 

Polytechnic  Institute  of  New  York 

Farmingdale,  NY  11735 

M.  H.  Bloom 

1 

Princeton  University 

Princeton,  NJ  08540 

F.  Hama 

1 

Purdue  University 

Lafayette,  IN  47907 

V.  W.  Goldschmidt 

1 

Rensselaer  Polytechnic  Inatitute 

Troy,  NY  12181 

R.  DlPrlma 

1 

L.  A.  Segel 

1 

Scripps  Inatitute  of  Oceanography* 

University  of  California 

La  Jolla*  CA  92037 

Dir. 

1 

Stanford  University 

Stanford*  CA  94305 

M.  van  Dyke 

1 

Stevens  Inatitute  of  Technology 

Hoboken,  NJ  07030 

J.  P.  Brealln 

1 

D.  Savitsky 

1 

Technion- Israel  Institute  of  Tech. 

Haifa*  Israel 

M.  Poreh 

1 

University  of  California 

Berkeley,  CA  94720 

Lib. 

1 

P.  Lieber 

1 

P.  Naghdi 

1 

W.  C.  Webster 

1 

J.  V.  Wehauaen 

1 

University  of  California*  San  Diego 

La  Jolla,  CA  2093 

J.  W.  Miles 

1 

University  of  Delaware 

Newark,  DE  19711 

J.  Wu 

1 

University  of  Hawaii 

Honolulu,  HI  96822 

J.  P.  Craven 

1 

University  of  Houston 

Houston,  TX  77004 

F.  Husaaln 

1 

University  of  Illinois 

Urbana,  IL  61803 

J.  M.  Robertaon 

1 

University  of  Iowa 

Iowa  City,  IA  52242 

J.  F.  Kennedy 

1 

L.  Landvtber 

1 

Univeraity  of  Maryland 

College  Park,  MU  20742 

J.  M.  Burgera 

1 

Prof.  Pai 

1 

University  of  Michigan 

Ann  Arbor,  MI  48105 

T.  F.  Ogilvie 

1 

C.  S.  Yih 

1 

University  of  Minnesota 

Minneapolis,  MN  55414 

Straub  Lib. 

1 

E.  Sllberman 

1 

University  of  Southern  California 

Loa  Angelea,  CA  90007 

J.  Laufer 

1 

Univeraity  of  Virginia 

Charlottesville,  VA  22903 

E.  J.  Gunter*  Jr. 

1 

Univeraity  of  Washington 

Seattle,  WA  98105 

A.  Hertzberg 

1 

C.  E.  Pearson 

1 

Virginia  State  College 

Petersburg,  VA  2 3803 

K.  M.  Agrawal 

1 

distribution  of  this  document  within  the  Applied  Physics  Laboratory  hat  been  made  in  accordance  with  a list  on  file  in  the  APL  Technical  Publications  Group. 


54 


INITIAL  DISTRIBUTION  EXTERNAL  TO  THE  APPLIED  PHYSICS  LABORATORY4 


ORGANIZATION 


LOCATION 


ATTENTION 


MISCELLANEOUS 

Australian  Embassy 
Engineering  Societies  Lib. 

Inatltuto  de  Tecnologls  Naval 
R.  H.  Kralchnsn 

Library  of  Congress.  Scl.  and  Tech.  Dlv. 
Maritime  Administration 

National  Bureau  of  Standards 
National  Research  Council 
National  Science  Foundation 
Society  of  Naval  Architects  and  Marine 
Engineers 

Webb  Institute  of  Naval  Architecture 


Washington,  DC  20036 
New  York,  NY  1001? 
Buenos  Aires,  Argentina 
Dublin,  NH  03444 
Washington , DC  20540 
Washington,  DC  20235 

Washington,  DC  20234 
Ottawa  7,  Canada 
Washington,  DC  20550 

New  York,  NY  10006 
Glen  Cove,  NY  11542 


Def.  Rea.  A Dev.  Attache 


R.  A.  Bastlanon 


Off.  Rea.  A Dev. 
Dlv.  Ship  Dea. 

G.  Kulln 

Lib. , Aero.  Lab. 

Eng.  Dlv. 


Tech.  Lib. 
E.  V.  Lewis 


'Initial  distribution  of  this  document  within  die  Applied  Physics  Laboratory  has  bean  mads  In  accordance  with  t list  on  fils  in  the  APL  Technical  Publications  Group. 


